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Abstract 

The generation of solar non-axisymmetric magnetic fields is studied based on a linear 
— n dynamo model in a rotating spherical frame. The model consists of a solar- 
like differential rotation, a magnetic diffusivity varied with depth, and three types of 
a-effects with different locations, i.e. the tachocline, the whole convective zone and 
the sub-surface. Some comparisons of the critical a-values of axisymmetric (m = 0) 
and longitude-dependent modes (m = 1,2,3) are presented to show the roles of 
the magnetic diffusivity in the problem of modes selection. With the changing of 
diffusivity intensity for the given solar differential rotation system, the dominant 
mode possibly changes likewise and the stronger the diffusivity is, the easier the 
non-axisymmetric modes are excited. The influence of the diffusivity and differential 
rotation on the configurations of the dominant modes are also presented. 
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1 Introduction 



Since Carrington (1863) suspected that sunspots did not form randomly over 
solar longitudes, many researches were undertaken to inv estigate the preferred - 



longi tudes of solar mag netic field. 'Active Longitu des' (ISvestka and Simon 



19691 ). 'Sunspot Nests' fICastenmiller et al. 1 . 119861 ). 'Hot Spots' flBai I . Il987f ) 



and so on were defined and presented. The well-k nown 'flip- flop' phenomena 



of lat e-type stars were also identified on the Sun (IBerdyugina and Usoskin 



20031 ). These more and more observational evidences indicate the involvement 
of non-axisymmetric, i.e. longitude-dependent large-scale magnetic fields in 
the formation and evolution of solar activities. 



Preprint submitted to Elsevier 



5 Eebruary 2008 



Although many papers reported the evidences for the existence of longitu- 
dinal inhomogeneities in sunspots distribution, the results were always in- 
consistent with each other on the number o f the a ctive longitudes, i.e. lon- 



gitudinal wave number m. lAltschuler et al. I (119741 ) analyzed the large-scale 



photospheric magnetic field in terms of surface harmonics a nd showed that 



the m odes m < 3 were prominent throughout the solar cycle. Ide Toma et al 



(I2OOOI ) argued that the total num bers of active regi o ns pe r rotation per hemi- 



sphere varied between and 7. ISong and Wang I (120051 ) suggested that 55 
percent of the solar magnetic flux can be represented by a mode of m = 6 
in the northern hemisphere and by m = 5 in the southern one. However, the 



identification of flip-flops on the Sun (IBerdyugina and Usoskin 1 . 12003| ) indi- 



cates the mainly non-axisymmetric mode m = 1. The ambiguous longitudinal 
wave numbers may arouse a series of questions. Which kinds of physical mech- 
anisms are responsible for the preferred mode to excite? For each mode, what 
is its configuration and which ingredients have the primary infiuence on the 
configuration? These are the main objectives of the paper. 

In order to explain the questions about the non-axisymmetric field of the Sun, 
it is natural to set up the non-axisymmetric solar dynamo m odel. There are 
two basic processes in the traditioii al a — Q dynamo model (jParker I . Il955l : 



Steenbeck. krause and Radler . 19661 ): (i) the generation of toroidal field by 
shearing a pre-existing poloidal field (the f2-effect) and (ii) the re-generation 
of poloidal field from the toroidal field (the a-effect) . Moreover, meridional 
circulation is always considered to transport the magnetic fiux. The a — Q 
dynamo is the limit of — Q dynamo in which the a-effect is included as 
a second generation source of the toroidal field. The role of the a-effect and 
the f2-effect in dynamo process has been well-researched and well-known. But 
there is another important ingredient, i.e. the magnetic diffusivity. Except for 
opposing the generation of magnetic field, what other roles does it play in the 
non-axisymmetric dynamo process? 

Magnetic diffusivity consists of molecular diffusivity and turbulent one. Since 
turbulent diffusiv i ty is m uch stronger, we neglect the molecular diffusivity. 
Yoshimura et al. I (Il984d ) showed that anisotropic turbulent diffusivity is a key 
factor on parity s election especially for th e latitudinal-gradient-dominated dif- 
ferential rotation. IChatterjee et el. I (120041 ) adopted different diffusivity-profiles 
(scalars) for toroidal and poloidal fields in a two-dimensional solar dynamo 
model, based on the helioseismically determined solar rotation profile and 
Babcock-Leighton-type a-effect. The solar-like dipolar parity was obtained. 
Obviously, the diffusivity has directly relation with the parity selection in dy- 
namo process. In the following, we also regard the diffusivity as a scalar and 
discuss its influence on the mode selection and magnetic conflguration during 
the non-axisymmetric dynamo process. 



In the next section, we set up a linear non-axisymmetric a'^ — Q dynamo model 
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in a rotating frame with solar inner core and introduce the spectral method 
to numerically solve it. We present the role of turbulent diffusivity as well as 
the differential rotation by the comparisons of the critical a-value of axisym- 
metric and longitude-dependent modes, and further show their configurations 
in Section 3. In Section 4, We give some discussion and conclusions. 



2 The Non-axisymmetric Dynamo Model 

2. 1 The basic equations and numerical scheme 

We solve the standard mean field dynamo equation 
dB 



V X [Ux B + aB -rjV X B] 

C/ o 



in the spherical polar coordinates (r, 6, (p). For the flow field U, only the (dif- 
ferential) rotation Q without the meridional circulation is considered for sim- 
plicity. 1] is the magnetic diffusivity. We non-dimensionalize the equation (1) 
in terms of the length Rq and the time Rq/tJo. Thus Eq. (1) becomes 

r) R 

— = RnVx{U'.xB) + i?„V x (a'B) - V x (r^'V x B), (2) 



where Rn 



OinRp) 



R, 



° — . rjo, ao and \^l\o are the reference values of 



Vo Vo 

the diffusivity in the convective zone (CZ), the a-effect and the differential 
rotation respectively. The quantities R^ and Rq are the dynamo numbers 
measuring the relative importance of inductive versus diffusive effects. 



At two interfaces r = Vo = 1.0 and r = = 0.6 (after non-dimensionalization) , 
both magnetic field and the tangential electric field must be continuous. The 
exterior r > 1.0 is a vacuum and eigensolutions are matched to a potential 
field. The radiative core is assumed to behav e as a perfect conductor (See 
Schubert and Zhang I (1200 ll ). Ijiang and Wang I (120061 ) for detail). 



Since the magnetic field is divergence-free, we expand B in term of two scalar 
functions h and g which represent the poloidal and toroidal potentials as 



flMoffatt I . Il978f ) 



B = V X V X rh{r, 9,(p,t) + V x rg{r, 9, (p, t). 



(3) 
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With given boundary conditions and linear frame, we can look for the eigen- 
solutions of the equations about h and g with the form 



[h{r,e,M,g{r,e,M] = [h{r,e,<P),g{r,e,<P)]e 



St 



(4) 



where s is the eigenvalue and can be written as s = a + iuj . Only the solutions 
that neither grow nor decay, i.e. the onset of dynamo actions {a = 0) are 
considered in the paper. For given Rq, the corresponding dynamo number Ra 
is the critical Ra and the corresponding value of ao is the critical a-value. 



Different modes are decoupled in the linear theory. We expand h and g at 
the onset of the dynamo action in terms of Chebyshev polynomial Tn{r) and 
surface harmonics pj^Q^^'t' in the meridional circular sector r G [0.6, 1.0], 9 G 
[0.0, tt] for given m as follows: 



N L 



im<f> 



n=0 l=m 
N L 



g = j:J2<iTn{ar - b)Pr{cose)e 



im<f> 



(5) 
(6) 



n=0 l=m 



where ar — b G [—1, +1]. and L are the truncation levels to get convergence. 
It varies with different dynamo numbers and different Q, r], a profiles, ; and 
; are eigenvectors. 



As pointed out by llvanova and Ruzmaikin I (Il985l ). the system of equations 



about h and g may be decomposed into two subsystems, i.e. odd or even parity 
with respect to the equatorial plane denoted by A and S. With the parame- 
ters adopted in our model, both odd and even parity solutions have nearly the 
same critical Ra and actions for given mode m. It should be enough to present 
one kind of parity for each mode. Since the odd parity for axisymmetric mode 
AO has been confirmed by the observations, we will discuss it and the oppo- 
site parity (even ones) for the non-axisymmetric modes (Sm). Note that, it 
does not mean there have no possibility for 5*0 a nd Am existing on the Sun . 
Through studying the heliospheric magnetic field, iMursula and Hiltulal (120041 ) 
gave the evidence fo r 5*0 mode on the Sun. For the non-axisymmetric modes, 
just hke IStix I (Il97ll ). only the modes m = 1,2,3 {SI, S2, S3) are considered 
since they are the most interesting ones in view of observational evidences. 
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Fig. 1. Left: Contours of constant angular velocity obtained from Eq. (7). The short, 
perpendicular ticks mark the 'downhill' direction. The dashed line denotes the center 
of the tachocline locating at 0.7Rq. The differential rotation is symmetric about the 
equator and only quadrant is shown. Right: Radial distributions of the three types 
of a-effects and the turbulent diffusivity rj. Values of y-axis are not scaled. 

2.2 Three ingredients in the model 



We use the following expression for the solar internal rotation (ISchou et al. 



19981 : ICharbonneau et al. I . Il999h 



2 w 



(7) 



where fls{d) = ^eq+o,2 cos^ 6'+a4 cos"^ 9 is the surface latitudinal rotation. The 
parametric values are set as Vc = 0.7Rq, w = O.OSi?©, i7c/2vr = 430.0 nHz, 
nEQ/271 = 455.8 nHz, a2/2n = -51.2 nHz, 04/271 = -84.0 nHz. It closely 
resembles the best-fit helioseismic solution. Left of Fig. 1 is the contours of 
constant angular velocity. Within the tachocline, rotation increases with dis- 
tance from the core at low latitudes, while it decreases at high latitudes. At 
intermediate latitudes (near 35°), rotation is almost independent on depth. 
Furthermore, we base our model on the rotating spherical systems with the 
rotation velocity Qc of the inner core. Thus the differential rotation in the 
rotating frame flc is 

n\r, 9) = (27125.8) x 1[1 + er/(2^^^)] 

2 w 

(1. - 1.98 cos^ 9 - 3.26 cos^ 9) (nHz). (8) 



The surface equatorial differential rotation is (27r25.8) nHz and we regard it 
as the characteristic value of the differential rotation \Q\o. Hence, Rq should 
8x1 0^*^ m'^s~^ 

be equal to ri — decided only by the reference value of the 
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diffusivity r], 



We make the simple assumptions that both a-effect and magnetic diffusivity 
T) are scalars. Whilst bearing in mind that they should be written as ten- 
sor quantities in a more realistic model. We use the analytical expression of 
Dikpati and Charbonneaul (119991 ) for the diffusivity profile as 



V{r) =Vc + ^[1 + erf {2- — 



)], 



(9) 



which can be seen in the right of Fig. 1 (solid line). The diffusivity rjo in CZ, 
is dominated by its turbulent contribution. In the stably stratified core, the 
diffusivity rjc is much lower because of the much less turbulence. In the fol- 
lowing, we take ?7c/''?o=0.01. The transition from high to low diffusivity occurs 
near the tachocline, which is coincident with the rotational shear layer. Here, 
rjo is far less definite but is widely known that it ranges from 2 x 10^'' cm^s~^ 
to 2 X 10^2 cm^s-^ 

The a-effect cannot yet be determined from observations. We consider three 
typical cases with different regions it works. They are in the tachocline, the 
convective zone and the near-surface. We denote them as CI, C2 and C3 
respectively. The expression for them can be written as 

a(r, e) = ao\[l + erfC-^)] ^[1 - erfC^)] cos9. (10) 



Ta, Tb and d change with different types of a-effects. ao occurred in the dynamo 
number Ra also presents the amplitude of the a-effect. The radial profiles 
are shown in the right of Fig. 1. The common angular dependence cos 6' is 
adopted, which is the simplest guaranteeing antisymmetry across the equator. 
Moreover, we do not consider the a-quenching since only the linear solutions 
are sought. 



3 Numerical results 



We perform some numerical explorations with the above three different a- 
profiles. In order to reveal the ingredients contributing to the mode selection, 
we test Rn at the values of 500, 2000 and 5000 and accordingly, the typical 
values of diffusivity are 1.6 x 10^^ cm^s"^, 4.0 x 10^^ cm^s~^, 1.6 x 10^^ cm^s~^. 
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Table 1 

The results of different modes for the case CI with different Rq and ?7o(10^^ cm^s~^). 
Middle three columns show critical a- values (cm s~^). The right three columns are 
the corresponding period (years) and phase shift between AO and Sm. Bold is used 
when the mode has the smallest a-value. 



Critical a Period T k. Phase shift 





500 


2000 


5000 


500 


2000 


5000 


Vo 


ib.U 


A n 
4.U 


i.D 


ib.U 


4.U 


i.D 


AO 


104.8 


9.91 


1.6 


1.85 


6.8 


16.88 


SI 


74.85 


6.70 


2.1 


12.18 


286.9 


7347.7 










0.847r 


O.GQvr 


1.047r 


S2 


73.60 


10.77 


4.2 


35.33 


1207.6 


3764.5 










1.307r 


0.167r 


1.337r 


S3 


82.56 


15.55 


6.4 


71.75 


7409.4 


1742.5 










0.487r 


0.797r 


1.707r 



Table 2 

The same as Tab. 1, but for Case 2. 



Critical a Period T & Phase shift 





500 


2000 


5000 


500 


2000 


5000 


Vo 


16.0 


4.0 


1.6 


16.0 


4.0 


1.6 


AO 


88.0 


5.89 


0.95 


4.55 


17.31 


43.06 


SI 


79.39 


8.50 


2.53 


26.06 


349.5 


2170.3 










0.727r 


0.737r 


0.747r 


S2 


83.47 


12.97 


4.65 


64.88 


739.2 


8024.5 










1.757r 


l.OOvr 


l.OOvr 


S3 


93.57 


17.7 


6.87 


119.6 


1524.6 


22587.4 










1.057r 


0.957r 


0.427r 



3.1 The comparisons of the critical a-values among the different modes 

The middle three columns of Table 1-3 show the critical a-values for different 

Rq, and r]o of the case CI that the a-effect concentrates in the tachocline 
(ra=0.675, rfe=0.725 and (i=0.025), C2 that the a-effect exists through the 
entire CZ (ra=0.7, rb=1.0 and (i=0.05) and C3 that the a-effect locates at the 
top of the surface (ra=0.95, r;,=1.0 and d=0.025), respectively. The modes 
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Table 3 

The same as Tab. 1, but for Case 3 and critical a-values (m s~^). 



Critical a Period T & Phase shift 





500 


2000 


5000 


500 


2000 


5000 


Vo 


16.0 


4.0 


1.6 


16.0 


4.0 


1.6 


AO 


26.17 


5.99 


1.63 


5.79 


1.16 


2.0 


SI 


16.77 


4.55 


1.88 


0.297 


0.293 


0.292 










O.lOvr 


I.IOtt 


0.77r 


S2 


17.84 


4.99 


2.11 


0.150 


0.147 


0.146 










O.lOvr 


O.IOtt 


1.297r 


S3 


19.12 


5.46 


2.22 


0.101 


0.100 


0.099 










0.727r 


0.257r 


1.277r 



with the lowest critical a-values are shown in bold. 



The general rule that the larger m is, the more difficult the mode to excite 
is only satisfied when Rq is large enough and the corresponding diffusivity 
is lower enough. For example, when Rq = 5000 (r/o = 1.6 x 10^^ cm^s^^) 
for all of the three cases and Rq = 2000 (r/o = 4 x 10^^ cm^s"^) for C2 and 
C3, the models will favor the axisymmetric mode AO. However, when the 
turbulent diffusivity is strong [rjo = 1.6 x 10^^ cm^s~^ for the three cases), 
the non-axisymmetric modes (m = 2 for CI, m = 1 for C2 and C3) have the 
lowest critical a-values and will be the preferred modes. Hence, the dynamo 
number Rq is the decisive parameter to select the dominant mode. According 
\Q\ R^, 

to Rq = — - — —, since the differential rotation is given, only the diffusivity in 

Vo 

CZ contributes to the result and hence should be the key factor to the mode 
selection. 



When Rq is large enough {Rq ^ Ra)-, the a-effect as the generation mecha- 
nism of toroidal field can be ignored and the — VL model is at the a — f2 
dynamo limit. At the a — Vl limit, the axisymmetric mode (m = 0) which is 
the smallest wave number and has the largest scale field is least depressed by 
turbulent diffusivity and will be the preferred mode. When rjo increases and 
the corresponding Rq decreases, the ratio of Rq to Ra also decreases. The role 
of a-effect in generation of the toroidal field increases and the a^ — Vt model de- 
viates the a — Vt dynamo limit gradually. The more it deviates from the limit, 
the easier the non-axisymmetri c modes are excited, which is in conformity 
with the asymptotic solution of iBassom et al. I (120051 ). 



The right three columns of Table 1-3 show the periods and phase shifts for 
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CI C2 C3 




0.0 0.5 1.0 0.0 0.5 1.0 0.0 0.5 1.0 



Fig. 2. Contours of the toroidal field B^f, of the axisymmetric mode AO in a meridional 
plane with i?Q=5000 for the three cases with different locations of a-effects. Solid 
(dashed) contours correspond to positive (negative) magnetic field components. The 
dot-dashed lines locate at 0.7 Rq marking the center of the tachocline. The radial 
distributions of the contour lines are consistent with that of the a-effects. But the 
latitudinal locations which are around 55° latitude do not change with the a-profiles. 



different modes with different Rq and rjo- For the axisymmetric mode AO 
of the cases CI and C2, it can be seen from the Table 1 and 2 that the 
period is nearly the inverse ratio of diffusivity and hence, is dominated by the 
diffusivity. The periods of the non-axisymmetric modes are much larger than 
that of AO. Therefore they are weakly oscill ating or rather steady coniparing 



with AO. This is consiste nt with observations f|Berdvug:ina and Usoskin I . I2003I ) 



and theory speculation (iFluri and Berdyugina . The co-existing of the 



two kinds of modes are used to explain active longitudes, flip-flops and o ther 
non-axisymmetric phenomena (IMoss I . l2004l : iFluri and Berdyugina 1 . 120041 ). In 
contrast, the oscillatory period of case C3 is much les s than the other two cases . 



This is consistent with the independ ent results of ( llvanova and Ruzmaikin 



19851 : iBigazzi and Ruzmaikin I . |200J). The periods depend on the diffusivity 



much less especially for the non-axisymmetric modes. C3 is characterized by an 
a-effect which has no overlap with the shear layer. It probably is the reason 
to cause the difference. The meridional circulation neglected in our models 
should play much more important roles in C3 to determine the cycle. 



Phase relation between the surface axisymmetric and non-axisymmetric poloidal 
field can be observed. Ruzmaikin et al. (2001) showed that there are about 
TT-shift between the modes AO and SI and the nearly s ame phase between 
5*1 and 5*2. The global phases (ISchlichenmaier and Stix I . Il995l ). based on the 



maximum values of B,. over the whole latitude ranges are computed. The phase 
shifts presented in Tab. 1-3 are the phase difference between AO and Sm. The 
results dependent on {rjo) and different cases. When Rq is large enough 
and diffusion is small, the phase shift can match the observation well. 
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3.2 The configurations of the magnetic fields for the three cases 



Fig. 2 shows the contours of the toroidal field i?^ of the axisymmetric mode 
AO in a meridional plane for the three cases when is 5000. The centers of 
the contour lines concentrate at the nearly same high latitudes (about 55°) 
where have the strong radial shear. But the radial locations of the three cases 
are different. They are around the tachocline, throughout the whole CZ and 
near the surface respectively, which are consistent with the radial profile of the 
a-effects. Hence, the strong radial shear of the differential rotation and the a- 
effect play a dominant role in the form of the configuration of the axisymmetric 
fields. 

Fig. 3 presents the distribution of of the distinguished non-axisymmetric 
modes for the three cases when Rq is 500. For CI and C2, the magnetic field 
is highly concentrated in the lower part of the tachocline, where the diffusivity 
is small (see Eq. (9)). But it is different for C3 with the sub-surface a-effect. 
Moreover, for all the three cases, the non-axisymmetric fields are localized 
around 35° latitude where has the weakest radial shear (see left of Fig. 1). 

The physical reason for the different localizations of the axisymm etric and 



non-a xisymmetric modes can be understood from MHD theory (iMoffatt 



19781 ). The differential rotation affects the two kinds of fields in different 
ways. It distorts the axisymmetric poloidal field perpendicular to the axis 
of rotation to form the axisymmetric toroidal field. By the way, this needs 
the non-zero diffusivity to assure the development of net toroidal flux. For the 
non-axisymmetric ingredients, they will expelled from the strong differential 
rotating high-diffusivity region. 



4 Conclusions and Discussion 



By the comparisons of critical a-values among the modes m < 3 for a wide 
range of turbulent diffusivity and dynamo number Rq accordingly, we obtain 
that dynamo number Rq is the decisive parameter for the modes selection. 
When Rq is large enough, the dominant axisymmetric mode would be fa- 
vored. With the decreasing of Rq, different non-axisymmetric modes will be 
preferred. Since Rq is the ratio of the differential rotation to magnetic diffusiv- 
ity, both of the two ingredients should contribute to the modes selection. The 
role of differential rotation and its physical mechanisms have been well-known 
that the strong value prefers the axisymmetric mode. For the deflnite differ- 
ential rotation given by helioseismology, the comparisons approve of the effect 
of the diffusivity in the mode selection and indicate that it is just contrary to 
the differential rotation. In other words, the strong turbulent diffusivity favors 
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CI 02 C3 




Fig. 3. Contours of the toroidal field i?^ of the non-axisymmetric modes in a merid- 
ional plane with Rq=500 for the three different cases. The line styles are the same 
as in Fig 2. For CI and C2, the field concentrates in the low diffusivity (below the 
dot-dashed line 0.7Rq) and weak differential rotation (about 35° latitude). For C3, 
it has the same latitudinal location but concentrates near the surface. 

the non-axisymmetric modes for the given rotation system of the Sun. On the 
other hand, the dominant mode given by the observation can be an important 
index to infer the inner turbulent diffusivity of the Sun or a star. What we 
have done is based on a linear non-axisymmetric a"^ — Q dynamo model in 
a rotating spherical frame which incorporates the solar-like differential rota- 
tion, the magnetic diffusivity varied with depth and three possible cases of the 
a-effects. 



On the actual Sun, it is rather than a single preferred mode. The whole physics 
during the dynamo process is not only the interaction between magnetic field 
and turbulent velocity field but also the interaction between the various modes. 
The latter may lead to a stationary field amplit ude and to a certain spectral 



distribution, even if the a- value is supercritical (jStix I . Il97ll ). T he nonlinear- 



ities, such as a non-axisymmetric d istribution of the a-ef fect (iMoss et al. 



2002 : Bigazzi and Ruzmaikin . 2004) and the a-quenching ( Moss . 19991 ) 



in- 



duce the different modes coupled together. The non-axisymmetric enhance- 
ment of the underlying magnetic field cau ses in the clustering of sunspots to 



form active longitudes ( IRuzmaikin I . l200ll ) 



Since the magnetic diffusivity is a diffusion of a vector field, the diffusivity pro- 
file should be expressed by a tensor. Yoshimura et al. (1984a ,1984b, 1984c) 
adopted the tensors for the diffusivity and demonstrated that it had close re- 
lations with parity selection. With a scalar expression, the growth rates for 
the odd and even parity modes are nearly the same (Seenbeck and Krause 
1969; Deinzer and Stix 1971) which also been demonstrated by our code. But 
in the paper, the scalar diffusivity is still used for two reasons. One is for the 
simplicity. We do not refer to the parity problem and only adopt the odd par- 
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ity for axisymmetric mode and even parity for the non-axisymmetric modes. 
The other is hmitation of the observations. The measurements of magnetic 
diffusivity inside of the Sun are not possible nowadays. Hence it is difficult to 
give its detailed distribution. 

The configurations of the toroidal fields for the axisymmetric and non-axisymmetric 
modes are different. The former prefers to locate at the strong radial shear 
region and the latter favors the weak radial shear and low diffusivity regions. 
These are caused by the different roles that the differential rotation and diffu- 
sivity played in the two kinds of modes. Furthermore, the meridional circula- 
tion is not considered in the work. It play s an important role in the axisymmet- 



ric mode (INandy and Choudhuri 1 . 12002| ). It carries the strong axisymmetric 



toroidal field produced at the high latitudes to the low ones and produces the 
active regions there with the magnetic buoyancy. We will discuss its infiuence 
on the non-axisymmetric field in the coming papers. 
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